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We analyze event-by-event fluctuations of the transverse momentum in relativistic heavy-ion col- 
lisions at y / sjviv = 200 GeV in the framework based on the fluctuating Glauber-model initial con- 
ditions, event-by-event (3 + 1) -dimensional viscous hydrodynamics, and statistical hadronization. 
We use the scaled fluctuation measure {ApTiApTj)/{{pT})- The identified "geometric" mechanism 
of generating the transverse-momentum fluctuations from the initial size fluctuations, transmitted 
to the final statistical-hadronization phase with hydrodynamics, is capable of easily reproducing 
the magnitude of the effect and explains the basic features of the data. On the other hand, it is 
somewhat too strong, hinting on modification of the popular Glauber approach to the earliest phase 
of the collision. We have checked that the considered measure is insensitive of the values of the 
shear and bulk viscosity coefficients, the freeze-out temperature, and the smoothing parameter for 
the initial distribution. It remains unaltered in the core-corona picture and is insensitive to the 
transverse-momentum conservation, approximately imposed in the statistical hadronization. 
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I. INTRODUCTION 

In Ref. [l| a new mechanism for generating the 
transverse-momentum fluctuations in relativistic heavy- 
ion collisions was identified. It is based on the ran- 
dom event-by-event fluctuations of the initial size of the 
formed system, its subsequent hydrodynamic evolution, 
and statistical hadronization. In the present work we fur- 
ther explore and extend this analysis, applying (3 + 1)- 
dimensional [(3 + 1)-D)] viscous event-by-event hydro- 
dynamics. The basic idea of Ref. [lj is as follows: Even 
when we consider a very narrow centrality class of events, 
e.g., with a strictly fixed number of wounded nucleons, 
N w , the size of the initial fireball fluctuates event- by- 
event due to the random nature of the nuclear collision 
in the Glauber treatment. These fluctuations are then 
transferred by hydrodynamics to the fluctuations of the 
generated transverse flow velocity. At freeze-out, this 
translates into the event- by- event fluctuations of the av- 
erage transverse momentum of hadrons produced in the 
event, (pr)- In essence, via simple scaling arguments, a 
more squeezed initial condition leads to more rapid ex- 
pansion, larger velocity flow, and higher (pt), while a 
swollen initial condition leads to slower expansion, lower 
flow, and lower (pt)- We will now explore this mech- 
anism through the use of state-of-the-art tools, such as 
GLISSANDO Q Monte Carlo code for the Glauber phase, 
(3 + 1)-D event- by- event viscous hydrodynamics [3|,|j] for 



the dynamical evolution, and THERMINATOR @, @| for the 
statistical hadronization at freeze-out. 

The event-by-event (pt) fluctuations in relativistic col- 
lisions have been actively studied theoretically [7-26] and 
experimentally |27H38| . as they may reveal relevant de- 
tails of the dynamics of the system, more accurate than 
contained in the one-body observables. Moreover, they 
are expected be sensitive to the critical phenomena at 
the phase transition, providing an important probe for 
these effects. 

Throughout the paper we use the notation 
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to indicate averaging in a given event, and averaging of 
the single-event averages over all events, respectively. 

The structure of the paper is as follows: in Sec. [TTJ we 
give the details of the Monte Carlo simulations of the ini- 
tial phase, focusing on the size fluctuations, Sec. UTT] pro- 
vides some necessary description of the applied (3 + 1)-D 
viscous hydrodynamics, while the statistical hadroniza- 
tion is described in Sec. lIVl We then proceed in SeclVlto 
presenting the results, which are compared to the data 
from the STAR and PHENIX collaborations. We inves- 
tigate the influence of model details on the results of our 
calculation, finding them very robust. In particular, the 
STAR measure of the event-by-event transverse momen- 
tum fluctuations is insensitive to the medium viscosity, 
freeze-out temperature, or the smoothing parameter of 
the initial distribution of sources. Our final conclusions 
and discussion is contained in Sec. I VII 



II. INITIAL STATE FLUCTUATIONS IN THE 
GLAUBER APPROACH 

The initial condition for hydrodynamics may be ob- 
tained from the Glauber approach, leading to the success- 
ful wounded- nucleon picture [39|,|40| (a wounded nucleon 
is a nucleon that collided inelastically at least once) or 
its descendants, such as the mixed model @, E|- When 
the initial condition is obtained via Glauber Monte Carlo 
simulations, the distribution of sources (wounded nucle- 
ons or positions of binary collisions) in the transverse 
plane fluctuates, reflecting the randomness in positions 
of the nucleons in the colliding nuclei. This leads to fluc- 
tuations of shape. 

The event-by-event fluctuations of the elliptic compo- 
nent of initial shape have been actively studied, as they 
lead to significantly enhanced elliptic flow |26|, I43 - I53J . 
They also generate odd Fourier components, absent from 
the event-averaged studies, such as the triangular defor- 
mation [54456J I , as well as higher-order components of the 
flow. Other interesting phenomena appear as the result 
of fluctuations, e.g., the torque effect [57| of the reaction 
planes at forward and backward pseudorapidities, or the 
directed flow at central rapidity [58|, [59| . 

We now describe in some detail the implementation 
of the Glauber model used in this work. The density of 
charged particles per unit of pseudorapidity, as a function 
of centrality, can be parametrized using a formula [4lll60l 
61] incorporating an admixture of binary collisions, iV^m, 
into the wounded-nuclcon model in the following way: 
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where a is a phenomenological parameter, a = 0.145 
for the highest RHIC energy of ^/snn = 200 GeV 
[6l|. The initial-state simulations are carried out with 
GLISSANDQ [ty, including a component from binary col- 
lisions. The parameter a in the initial distribution is 
somewhat smaller from the value extracted from the fi- 
nal distributions (see the following). The difference is due 
to the longitudinal expansion and entropy production in 
the (3 + 1)-D viscous hydrodynamic expansion [jj. 

The positions of nucleons in each of the colliding nu- 
clei are randomly generated from a Woods-Saxon dis- 
tribution, with an additional constraint enforcing the 
short-range repulsion, namely, that the centers of nu- 
cleons in each nucleus cannot be generated closer than 
the expulsion distance d = 0.9 fm. Nucleons from the 
two colliding nuclei are wounded, or a binary collision 
occurs, when their centers get closer to each other than 

the distance ycr^^/n, with a^iv denoting the inelas- 
tic nucleon-nucleon cross section. For the highest RHIC 
energy of ^/sWn = 200 GeV one has a^% = 42 mb0 



1 One may more appropriately use a Gaussian wounding profile in- 
stead of the applied hard-sphere wounding profile, but the results 
do not differ significantly in the case of size fluctuations 1621 . 



The notion of sources, originally limited to the trans- 
verse plane, may be extended on the rapidity depen- 
dence of the particle emission. Although this extension 
is not crucial for the present study, focused on the mid- 
rapidity region, we include it for the integrity of the pa- 
per. The spatial pseudorapidity fa||) distribution of the 
emission profile is given as the sum of contributions from 
the forward- and backward-moving wounded nucleons. 
Within such an extended framework Bialas and Czyz 
have properly described [63J the pseudorapidity distribu- 
tions of charged particles in the d— Au collisions. There- 
fore, we assume an asymmetric emission profile [63L I64J 
peaked in the forward (backward) rapidity for the for- 
ward (backward) moving wounded nucleons, denoted as 
/+fa||) (/-fa||)), 
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where j/bcam is the beam rapidity. The initial profile in 
space-time rapidity is 

/fan) = cx P (- {r,ll ~J o)2 e(\ m \ - m )\ , (4) 

with rjQ — 1.5, a , ; — 1.4 [3j. The initial entropy density 
is assumed to have a factorized form 

s(x,y, m ) = Kj2f±(m)9i(x,y) [(l-a)+N9° u a] . 

i 

(5) 

Here Nf° n is the number of collisions of the participant 
nucleon i, and 
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implements a Gaussian smearing, replacing the point-like 
source at the transverse position (xj,j/i) with a Gaus- 
sian profile. The smearing parameter is taken to be 
w = 0.4 fm, and the overall scale factor is n — 2.5 GeV. 
The parameter a of the mixed model is fixed to reproduce 
the dependence of dN/dij on centrality. In the (3 + 1)- 
D viscous hydrodynamic model the optimum value is 
a = 0.125 at the top RHIC energies Q. We remark that 
the mixed model works also very well for the description 
of multiplicities a the LHC energy of y / J/v]v = 2.76 TeV, 
where a = 0.15 [H, HI]. 

In Fig. Q] we show two snapshots of typical config- 
urations of sources in the transverse plane generated 
with GLISSAND0. The dots indicate the positions of the 
wounded nucleons. Since we have in mind the distri- 
butions as starting conditions for the event-by-event hy- 
drodynamics, we need to smear out the point-like dis- 
tributions. The smearing procedure, although physically 
motivated and necessary, is somewhat arbitrary in intro- 
ducing a smearing scale. In Fig.[TJ the contours show the 
smeared entropy density, s, with w = 0.4 fm. Although 
both selected events correspond to the same number of 



6 


N„= 


=100 ^-s 


<r 2 > 1/2 =3.14fm 




/T* 


^\ <p T >=563MeV 






\W 




£ 3 








*" 




/^Dr 


f&*) 










-3 


a) 






-6 




(•/ 






3 
x [fm] 

FIG. 1. (Color online) Two typical configuration of wounded 
nucleons in the transverse plane (dots) generated with 
GLISSANDO and the corresponding contours of the smeared 
density of entropy, s. Solid, dashed, and dotted lines cor- 
respond to isentropes at s = 0.05, 0.2, and 0.4 GeV -3 , re- 
spectively. The densities for the two events have radically 
different r.m.s. radii of 3.14 and 2.38 fm, respectively, despite 
the equal number of the wounded nucleons, N w = 100. 



wounded nucleons, N w = 100, they have radically differ- 
ent r.m.s. radii, which after the hydrodynamic expansion 
results in different transverse flows. 

To have a simple size measure we look at the average 
transverse size of the initial fireball, defined in each event 
via the mean squared radius at the central space-time 
rapidity 



(r 2 )^ 



_ J dxdy(x 2 + y 2 )s(x, y, 0) 



J dxdy s(x, y, 0) 
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In the following we use the notation (r) = (r 2 ) 1 / 2 . The 
point, clearly seen from Fig. [1] is that even at precisely 
fixed centrality the size (r) fluctuates [lj. The feature 
is presented quantitatively in Fig. [3J where we plot the 
event-by event scaled standard deviation of (r) obtained 




FIG. 2. (Color online) Event-by-event scaled standard devi- 
ation of the size parameter (r), evaluated at fixed values of 
the number of wounded nucleons N m from the initial entropy 
density for events used in hydrodynamic simulations. 



at each N w . As expected, c((r))/((r)) is a decreasing 
function of N w . 

As noted in Ref. [l|, very similar curves to Fig. [2] arc 
obtained for other variants of Glauber models, such as 
models with overlaid distributions of particles produced 
from the sources J2|, simulations applying a Gaussian 
wounding profile [67| for the A^A~ collisions, or the use 
of the nucleon distributions including realistic (central) 
AW correlations of Ref. [68-70]. This means that the be- 
havior of the initial geometry shown in Fig. [5] is robust, 
essentially reflecting the statistical feature of the Glauber 
approach. 



III. VISCOUS EVENT-BY-EVENT 
HYDRODYNAMICS 

It is widely believed that a successful and uniform de- 
scription of the physics of relativistic heavy-ion collisions 
is achieved with the help of relativistic hydrodynamics 
(for reviews see, e.g., [7ll473l ]). Event-by-event hydrody- 
namic calculations for fluctuating initial conditions have 
been performed for perfect fluid [U [H, [H, [zHzl] and 
for the viscous case [4J, |77H79| , focusing on collective flow. 

In the second-order viscous hydrodynamic formalism 
the hydrodynamic equations 
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with the energy-momentum tensor 

T^ v = (e + p)u"u u - vf v + Tr"" + nA"" (9) 

are supplemented with equations for the stress correc- 
tions from the shear, 
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and the bulk viscosity, 
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Here V M = A^c*,,, while r\ and £ denote the shear and 
bulk viscosity coefficients, respectively. In our default 
calculations we use constant r\j s — 0.08, (,/s — 0.04 in 
the hadronic phase, r n ~ 3i]/(Ts), and rn = 7V. To 
test the sensitivity of our results on viscosity, we perform 
calculations for -q/s = 0.16, £/ ' s — 0.04 and -q/s = 0.08, 
C/s = 0.08 as well. 

The applied equation of state is a crossover equation 
of state, interpolating between the lattice-QCD results 
at high temperatures [83j and a hadronic gas equation 
of state at low temperatures. The construction of the 
equation of state follows the method of Chojnacki and 
Florkowski [H[ (for details see [!]). 

In this work we apply the event-by-event (3 + 1)- 
D viscous hydrodynamics [J, l77|, starting the evolu- 
tion at 0.6 fm/c. The configurations of wounded nu- 
cleons and binary collisions corresponding to the central- 
ity range - 70% are generated with GLISSANDO. The 
procedure does not fix the impact parameter for each 
centrality bin, as the Monte-Carlo scheme picks the im- 
pact parameter in each event according to the distribu- 
tion P(b) = da ine i(b)/(dba ine i) 085]. For each configu- 
ration of wounded nucleons a hydrodynamic evolution is 
calculated starting from the density ([5]). 
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FIG. 3. (Color online) Averaged inclusive transverse momen- 
tum vs. number of wounded nucleons, N w . The data (ex- 
trapolated to the whole pr range) come from the PHENIX 
Collaboration [92|] and show the charged pions (down trian- 
gle), charged kaons (up triangles), and protons and antipro- 
tons (squares) . The lines correspond to our model calculation 
with r)/s = 0.08, C/s = 0.04, and T f = 150 MeV (solid lines), 
and T f = 140 MeV (dashed lines). 



respectively, with /q denoting the equilibrium distribu- 
tions and c s standing for the velocity of sound. In the 
local rest frame the normalization constant is 



IV. STATISTICAL HADRONIZATION 

The last stage of our approach is the simulation of the 
statistical hadronization at freeze-out [861 ] (for a review, 
see, e.g. (zl]) with THERMINATOR [111]. The code includes 
all resonances and decay channels from SHARE 87]. The 
particles (stable and unstable, which subsequently de- 
cay) are formed at the freeze-out hypersurface according 
to the Frye-Cooper formula. In the case of viscous hy- 
drodynamics, the momentum distributions at freeze-out 
are modified by the viscous corrections. The shear and 
bulk viscosity corrections are [88[ 
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This is simply achieved by generating a uniform distribution in 
b 2 and accepting those events where at least one NN interaction 
occurred. 
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where the sum runs over all the hadron species. The 
(single-fluid) hydrodynamic evolution uses an equation of 
state with zero chemical potentials. However, the chem- 
ical potentials are reintroduced in the Frye-Cooper for- 
mula with the ratio \ijT fixed through the fits to the 
particle ratios at the chemical freeze-out, which works 
properly at the RHIC energies J91j . 

Before showing the pr-correlation results, let us 
state that our approach properly describes the relevant 
one-body features of the collisions, in particular, the 
transverse-momentum spectra. As an example, in Fig. [3] 
we show the inclusive average transverse momentum as 
the function of N w for pions, kaons, and protons and 
antiprotons for our default parameters Tf = 150 MeV, 
rj/s = 0.08, C/s = 0.04 (solid lines). The result compares 
favorably to the PHENIX data 92]. The agreement is 
important, as it shows that we have the correct one-body 
background to study correlations. Fixing the freeze-out 
temperature of Tf = 150 MeV reproduces the transverse 
momenta of identified particles at midrapidity. To check 
the sensitivity of the results of the freeze-out tempera- 
ture and viscosity, we have investigated also the cases 
when one of the parameter is modified from the default 



value to Tf = 140 MeV, rj/s = 0.16, or Q/s = 0.08. The 
calculations with a lower freeze-out temperature or with 
an increased shear or bulk viscosity give average trans- 
verse momenta within the range of the systematic er- 
rors quoted by the PHENIX Collaboration. Admittedly 
there is some model dependence on parameters, but it is 
weak, and the default parameters serve as an optimum 
choice. 

It has been noted that event-by-event hydrodynam- 
ics with lumpy initial conditions yields harder spectra 
than hydrodynamics starting with averaged initial con- 
ditions [93|. This effect follows from higher gradients 
in the lumpy initial condition. To compensate, i.e., to 
soften the spectra, one needs to run hydrodynamics for a 
shorter time, i.e., to higher freeze-out temperatures [3[. 



V. RESULTS 

The simulations presented in this section employ the 
experimental cuts in the STAR [32j (0.15 GeV < pr < 
2 GeV) and PHENIX J29J (0.2 GeV < p T < 2 GeV) anal- 
yses. In both cases \r/\ < 1. Our samples have 100 events 
at each considered centrality bin. These, involving the 
hydrodynamic evolution, are time-consuming to gener- 
ate. To increase the accuracy of the statistical hadroniza- 
tion, we generate 200 THERMINAT0R events for each hydro 
event. 

Our determination of centrality matches closely the ex- 
periment. In the case of STAR [32J, the multiplicity of 
generated charged particles in the window |r;| < 0.5 is 
used to determine the centrality bins. In the case of 
PHENIX [29J, where a combination of signals from the 
BBC and ZDC detectors is used, we simply take the num- 
ber of wounded nucleons N w as the variable fixing the 
centrality. 



A. Fixed number of wounded nucleons 

For better understanding, we begin the analysis for 
the event-by-event fluctuations by selecting a very nar- 
row centrality class, with N w = 100. We run GLISSAND0 
to generate the initial conditions, carry out our event- 
by-event hydrodynamics, and, finally, run THERMINATOR 
and compute (pt) in each event. As argued before [l|, 
the fluctuations of the initial condition manifest them- 
selves in the fluctuations of the initial size (r) . In Fig. @] 
we plot the values of (j>t), histogrammed in bins of (r). 
Each point corresponds to one event, while the bars give 
the event- by- event average, ((pr))- We note a clear an- 
ticorrelation of ((pt)) and ((r)). This shows that in a 
full-fledged event-by-event simulation the basis qualita- 
tive argument holds: for a squeezed initial the system 
expands with the larger flow velocity hand acquires a 
higher average transverse momentum, ((px)), than for 
the stretched state. The same effect can be observed 
when comparing case by case events generated from dif- 
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FIG. 4. (Color online) Averaged transverse momentum as the 
function of the initial size (r) for events with a fixed number 
of wounded nucleons, N w — 100. Viscous (3 + 1)-D event-by- 
event hydrodynamics with rj/s — 0.08, £/s = 0.04, and Tf = 
150 MeV is used. The scattered small dots show (pr) obtained 
in individual events, while the bars show the event-by-event 
averages {{pr}} in the selected bins of (r). The anticorrelation 
is apparent, with lower size (r) resulting in higher {{pt})- 



ferent initial conditions (Fig. [TJ. The event with a 
squeezed initial density has a larger transverse flow and 

(pt)- 

The fit to the histogram bars in Fig.@]yields ((pr)) = 
0.79 - 0.07((r)) GeV/fm, which in turn gives 
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in the considered range. This result can be written as 
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which may be compared to the estimate of Ref. [94[ , 
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with P and e denoting the average pressure and en- 
ergy density during the evolution of the system. Thus 
P/e ~ 0.15, which is the right ball park for the applied 
equation of state [95[ . 



B. Transverse momentum fluctuations vs. 
centrality 

Now we come to the main results of this paper. In 
order to compare to the data, we analyze the STAR cor- 
relation measure [32j, (ApxiA-PTj): defined as 
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FIG. 5. (Color online) Comparison of the theoretical predic- 
tions for (Ap Ti Ap T3 ) 1/2 / ((p T )) (for yfs^ = 200 GeV) to 
the experimental data extracted from the STAR Collabora- 
tion [33 (squares). The dots correspond to simulation with 
event-by-event (3 + 1)— D viscous hydrodynamics with our de- 
fault parameters T f = 150 MeV, rj/s = 0.08, (/s = 0.04. The 
statistical errors of the model simulation are obtained with 
the jackknife method. The experimental statistical errors are 
negligible. 
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FIG. 6. (Color online) Comparison of the theoretical pre- 
dictions for (Ap Ti Ap Tj ) 1/2 /((PT)) (for V«aw = 200 GeV) 
to the experimental data from the PHENIX Collaboration 
:29J (squares). The dots correspond to simulation with event- 
by-event (3 + 1)-D viscous hydrodynamics with our default 
parameters T> = 150 MeV, r)/s = 0.08, (/s = 0.04. The 
crosses indicate the approximate result from Ref. [lj for per- 
fect (2 + 1)-D hydrodynamics with averaged initial conditions 
from the mixed model. The statistical errors for the model 
simulations are obtained with the jackknife method. 



where N cv is the number of events, N k the multiplicity 
in event k, and 

N k N k 
C* = E £ (p*-«PT»)fe-((pT»), (19) 
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Introducing the mean momentum in event k, denoted by 
(pr)k, we can transform 
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(Ap Ti Ap Tj ) = — ^ var((pr)) 
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Thus the STAR correlation measure is the difference of 
two terms: one involving the variance of the mean mo- 
menta in events, and the other being the event-averaged 
variance of the momentum in each event decided by the 
multiplicity of this event. Note that expression (|22|) in- 
volves only single sums in a give n event. As a mat- 
ter of fact, the STAR analysis 32j replaces ((pt)) with 
the quantity (j»t) (-^charged), the average momentum as a 
function of the number of charged particles in the pseu- 
dorapidity bin |t/| < 0.5 - the same as used to determine 
centrality. The function is obtained by a numerical fit to 
the results prior to the analysis of the correlations. The 
method slightly reduces the value of (ApTiApxj) ■ We 
follow the same prescription. 

Our results are shown in Fig. [SJ where we compare 
the theoretical points (circles) to the experimental data 
from the STAR Collaboration [32[ (squares) . At low cen- 
tralizes, the model calculations overshoot the data by 
about 50%, yielding more pr fluctuations than needed. 
This conclusion supports the original findings of Ref. [l[ 
in the present state-of-the-art event-by-event treatment. 
We have checked for a few centrality bins that modify- 
ing the shear or bulk viscosity coefficients, the freeze-out 
temperature or the width of the smearing Gaussian in 
the initial conditions docs not change the results at the 
level of the statistical errors of our calculations (see sect. 

Ecj. 

Nevertheless, we note a proper magnitude of the effect 
and the correct dependence on centrality, Also, since the 
results of Fig. [5] very weakly depend on (Jnn fl], with 
the expectation that the hydrodynamic "push" is similar 
at different collision energies, our results should weakly 
depend on the incident energy. This is a desired feature, 
as the STAR data [32] are very similar from ^Jsnn = 
20 GeV to 200 GeV. 

The statistical errors of the model simulations in Fig. [5] 
are estimated with the jackknife method. Essentially, the 
relative error is equal to l/\/2n, where n = 100 is the 
number of the hydrodynamic events in the considered 
centrality class. 

The PHENIX Collaboration [H] published results on 
the ratio of the pr fluctuations using the measure 
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where 
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and w™ is the same quantity obtained with mixed 
events. For small dynamical fluctuations and sharp dis- 
tributions in the multiplicity variable one can estimate^ 



(ApTiApTj) ~ 2F P 



var(p T ) 
' (N) : 



(25) 



where var(px) denotes the inclusive variance of the trans- 
verse momentum distribution, and (N) is the average 
multiplicity of the detected particles in the considered 
centrality class. The values of the quantities on the right- 
hand side of Eq. (J25J) are available from the PHENIX Col- 
laboration web page associated with Ref. [29f . We stress 
that the result Eq. (|2"5) is approximate, but sufficiently 
accurate [96] for our purpose. A more direct comparison 
to the PHENIX data could be achieved with the mixing 
technique, however, this is beyond our reach due to a 
very limited number of the model events. 

The result of the analysis is shown in Fig. [SJ with sim- 
ilar conclusions as from Fig. [SI i.e., the model points are 
above the experiment. We also show that the results 
of applying the event-by-event viscous hydrodynamics 
(dots) are very close to the approximate calculation of 
Ref. [l( for perfect (2 + 1)-D hydrodynamics with aver- 
aged initial conditions (crosses). 

The dependence of the fluctuation measure 
(ApTiApTj) 1 ' 2 /((pt)) on the upper transverse- 
momentum cut-off has been measured by the PHENIX 
Collaboration [29]. As can be seen in Fig. [JJ the 
fluctuations in the model increase with the cut-off, 
following closely the trend observed in the data. This 
cross-checks that the observed pr-dependence of the 
transverse momentum fluctuations can be interpreted as 
a hydrodynamic flow effect. 



C. Dependence on model parameters 

In this section we investigate the dependence of our 
predictions on the model parameters, such as the viscos- 
ity coefficients of the medium, the freeze-out tempera- 
ture, Tf, or the smoothing parameter, w. For this pur- 
pose we have run simulations with various values of these 
parameters at fixed N w = 100. 

In Fig. [8] we compare the dependence of (pr) on (r) for 
several variants of viscous hydrodynamics. The points 
correspond to mean ((pt)) in a given (r) bin, and the 



3 The relations between various popular correlations measures in 
this limit are discussed in the Appendix of Ref. [96| . One of 
us (WB) thanks Jeff T. Mitchell for the discussion concerning 
Eq. J23} 



A 
A 

V 
V 



A_ 
P 

< 

P 
< 

V 



3.5 ; 

3.0 ; 

2.5 ; 

2.0 ; 

1.5 : 

1.0 : 



[] 



II 



Pt** [GeV] 

FIG. 7. (Color online) Dependence of {Ap T iAp T j} 1/2 / ({pr}} 
(for ^/sjvjv = 200 GeV) on the upper transverse-momentum 
cut-off, compared to the experimental data from the PHENIX 
Collaboration [2SJ] for centrality 20 — 25% (squares). The dots 
correspond to simulation with event-by-event (3+l)-D viscous 
hydrodynamics with our default parameters Tf = 150 MeV, 
r\j s — 0.08, (/s = 0.04. The statistical errors of the model 
simulation are obtained with the jackknife method. The 
shaded band in the upper right corner represents the error 
band of the result of the simulation with an infinite upper 
momentum cut-off. 



curves are linear fits to these points. Our default re- 
sult is for the set of parameters rj/s — 0.08, (/s = 0.04, 
Tf — 150 MeV, and w = 0.4 fm, indicated with the up- 
triangles and dot-dashed line in the plot. We then do our 
comparison by changing one of the parameters: the shear 
or bulk viscosity coefficient, the freeze-out temperature, 
or the smoothing parameter. As expected, increasing the 
shear viscosity or decreasing the freeze-out temperature 
leads to a hardening of the spectra, i.e, higher ((pt))- On 
the other hand, increasing bulk viscosity leads to a reduc- 
tion of the effective pressure and a decrease of the average 
transverse momentum. Increasing the smearing width of 
the initial density distribution yields smaller gradients, 
and reduces the transverse push. The described behav- 
ior holds bin-by-bin in the (r) variable giving the size of 
the initial geometry. 

At the same time we note that the slope of the depen- 
dence of ((pt)) on (r) changes as well. For all the studied 
cases it turns out that to a high accuracy 



d(p 7 



d(r) (p 2 



0.31. 



(26) 



As a result, according to the arguments of Sec. IV Al 
the scaled measure (ApnApTj) / ((pt)) 2 is hardly modi- 
fied. Increasing the transverse pressure or the local gra- 
dients gives a larger transverse flow and, simultaneously, 
larger fluctuations, such that the scaled fluctuations of 
the transverse momentum practically do not depend on 
viscosity, the freeze-out temperature, or the smearing 
of the initial conditions. We thus find that the scaled 
PT-Auctuations are dominated by the fluctuations of the 
transverse size of the initial fireball. This feature makes 
the scaled measure particularly suitable for constraining 



M=100 



average transverse momentum, 
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(( PT )) = (l- C )((^))+c((^ OT0 )), 



(28) 



FIG. 8. (Color online) Dependence of {{pr)) on (r) for the 
variants of the hydrodynamic evolution: calculation with 
r)/s = 0.08, C/« = 0.04, T f = 150 MeV, and w = 0.4 fm 
(default) (triangles and dashed-dotted line), with rj/s = 0.16 
(filled circles and dotted line), with Tf — 140 MeV (squares 
and solid line), with f/s = 0.08 (open circles and dashed line), 
and with the width of the smearing w = 0.6 fm (stars and 
long dashed line). The points with the errors bars represent 
the histogram of the average momentum {{pr}) as function of 
the r.m.s radius of the initial density of the event. The lines 
represent linear fits to the points. 



the models of the initial phase. In a similar way, the di- 
rected flow at central rapidity has been proposed as a tool 
to limit the dipole deformation of the fireball predicted 
by different models of the initial state [97| . 



D. Other effects 

In peripheral collisions the particle emission can take 
place in the thermalized, collectively expanding core, as 
wel l as in the outer corona, where rescattering is small 
98 101]. In the following we estimate the transverse mo- 
mentum fluctuations in the case where the particles are 
emitted from these two sources, the core and the corona. 
Different definitions of the dense core are possible; we 
use the prescription that a wounded nucleon belongs to 
the core if it collides more than once. This choice of the 
separation between the core and the corona describes the 
centrality dependence of the strangeness productio n and 
of the effective slopes of the particle spectra [99|, Il02j j . 
The particle density at central rapidity is a sum of the 
contributions from the corona and the core, 



dN, 



charged 



dt] 



dN NN 

Lil v charged 

dr/ 



(N c , 



■PN corc ), (27) 



with iV C ore + -^corona = N w . To reproduce the central- 
ity dependence observed experimentally we choose the 
parameter (3 = 1.75, which effectively describ es th e en- 
hanced production in the thermalized matter [lOOj. The 



is a combination of the average momentum in a NN- 
collision and of the average transverse momentum of par- 
ticles emitted from the core. The fraction of particles 
emitted from the core is c = (3N cole /(f3N COTe + N COIona ,). 
By neglecting the fluctuations of the number of nucleons 
in the core, we get 



(Ap Ti Ap T j] 



-(l-c) 2 (Ap Tl A PT ^' 

+ <-<i-c) 2 («p£"»-«pS- 



X-'* corona/ 
„2/i „\2/ 



„core\\ \2 



e »)^ 



+ c 2 (Ap Tl Ap Tj )Y 



(29) 



The first term is a contribution from the N corona indepen- 
dent NN sources, the second term comes from the dif- 
ference of the transverse momenta from the two sources, 
and the third term is a contribution from the hydrody- 
namically expanding core. We use the PHENIX data 
[29( on the average transverse momentum and its fluctu- 
ations in the pp-collisions, and take for (ApTiApTj) colc 
the results from the hydrodynamic model calculation. 
With all elements of Eq. (|2"9")l combined, we find that 
(ApTiApTj} 1 ^ 2 1 ((pt)) is changed very little compared 
to the results of Sec. IV Bl For the most central collisions, 
where the corona contribution is tiny, naturally the effect 
is negligible. In peripheral collisions (c = 60 — 70%) all 
terms of Eq. (f2"9")l contribute to the transverse momentum 
fluctuations. However, in the last term the reduction of 
the core due to the (1 — c) 2 factor is compensated with in- 
creased fluctuations <r((r))/(r), such that with all terms 
combined the change is at the level of 10%. At interme- 
diate centralities the reduction effect is at a similar level. 
Therefore the core-corona model does not improve nor 
deteriorate the agreement with the experimental data. 

We have also checked that imposing a finite detector 
acceptance, by simply accepting a simulated particle with 
the typical probability of 50%, does not alter the results 
for (ApTiApTj) I ((pt)) ■ This is a feature of the scaled 
Pt fluctuation measure [29|, |37| • 

Finally, we have estimated the possible effect of the 
global transverse-momentum conservation, not imple- 
mented in the standard simulations of the statistical 
hadronization with THERMINAT0R. This can be approx- 
imately achieved by accepting only those events which 
have limited total transverse momentum. Specifically, 
we consider the quantity P 2 = (^2pi, x ) 2 + {J2Pi,v) 2 m a 
given event and include the event for the further analysis 
when P is less than a specified value, which is gradually 
decreased. With our statistics we are able to reduce the 
limit for P down to 15 GeV (there are a few hundred 
of particles in the event), which leaves about 5% from 
all (unconstrained) events for the most central case. No 
noticeable effect is detected, therefore the considered cor- 
relation measure is not sensitive to the global transverse- 
momentum conservation. 



VI. CONCLUSIONS 

The initial shape and the volume of the fireball fluctu- 
ate due to the random nature of the Glauber approach. 
As is well known, the subsequent hydrodynamic evolu- 
tion carries over the asymmetry of the shape of the fire- 
ball into anisotropies of the particle spectra. A similar 
mechanism, analyzed in detail in this work, transmits 
the event-by-event fluctuations of the transverse size of 
the fireball into the fluctuations of the average transverse 
momentum in each event, as identified in [l|. Here are 
the main findings of our analysis: 

1. The state-of-the-art event-by-event viscous (3 + 1)- 
D hydrodynamic calculations with fluctuating ini- 
tial conditions confirm that fluctuations of the 
mean transverse momentum in each event are gen- 
erated from the fluctuations of the initial geometry. 

2. The amount of scaled transverse momentum fluc- 
tuations is determined by the scaled fluctuations 
of the transverse size of the fireball. We observe 
an anticorrelation of the initial size of the fireball 
and of the transverse momentum generated in an 
event. The expansion of a source of larger extent 
yields smaller px then in the case of a squeezed 
source, and vice versa. 

3. Hydrodynamic expansion is applied to an ensemble 
of events, corresponding to centralities — 70% in 
Au-Au collisions at y/sjyN — 200 GeV. We find a 
similar magnitude and centrality dependence of the 
scaled momentum fluctuations (ApxiApTj) / ((pt)) 
as in the STAR [H and PHENIX experiments [29J| . 

4. The dependence of the results on the upper cut-off 
for the transverse momentum of the particles agrees 
nicel y w ith the data from the PHENIX Collabora- 
tion [29j . 

5. However, the initial density from the mixed model 
(wounded nucleons with an admixture of binary 
collisions), tuned to reproduce the particle multi- 
plicities, yields a visible overprediction of the ob- 
served value of (ApxiApTj) / ((pt}) in the whole 
centrality range. For most central events the over- 
prediction is at the level of 50%, while it gets rela- 
tively closer to the data with increased centrality. 



6. Hydrodynamic expansion yields a stronger trans- 
verse push and, simultaneously, stronger px- 
fluctuations when the shear viscosity is increased, 
the freeze-out temperature is lowered, or if the bulk 
viscosity is lowered. However, the scaled fluctua- 
tion measure (ApTiApTj) / {(pr)) shows very little 
changes with these modifications of the physical pa- 
rameters. 

7. Predictions of our approach remain essentially un- 
changed when the core-corona mechanism of par- 
ticle emission is incorporated. Other effects, such 
as the transverse-momentum conservation or the fi- 
nite detector acceptance do not affect the results, 
either. 



The above points indicate that the identified "geomet- 
ric" mechanism of generating the transverse-momentum 
fluctuations from the initial Glauber-like model is, on the 
one hand, very important, easily reproducing the size of 
the effect and catching the basic features of the data, on 
the other hand, it is somewhat too strong. That hints on 
an improvement of the popular Glauber approach of the 
initial phase. We recall that the calculations using the 
averaged initial conditions [l| show that the fluctuations 
of the initial size are reduced if the density in the fireball 
is determined with the wounded nucleons only, i.e., with- 
out the admixture of binary collisions. At the same time, 
however, the model with the wounded nucleons only fails 
to generate the proper multiplicity dependence on cen- 
trality, thus is less realistic. Moreover, hydrodynamic 
fluctuations in th e evolutio n could add another source of 
(pt) fluctuations [l03lll04| . Thus, it remains a challenge 
to understand in detail the earliest phase of the colli- 
sion and reproduce in a uniform way the rich collection 
of the one-body and the correlation data, including also 
the harmonic flow. 
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